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Abstract: We developed tetranucleotide-repeat microsatellite markers for the masu salmon 
(Oncorhynchus masou) complex. 454 pyrosequencing was used to discover repeat motifs, 
and seven polymorphic microsatellite-primer sets were identified. The number of alleles 
detected at each locus ranged from four to 24 and the expected heterozygosity varied from 
0.57 to 0.92. Cross-subspecies amplification for O. m. masou, O. m. ishikawae and 
O. m. subsp. was successful. These microsatellites can be utilized in studies of genetic 
structure, genetic diversity, and intra- and inter-subspecific hybridization, making a 
contribution to conservation and management of the Oncorhynchus masou complex. 
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1. Introduction 

The masu salmon (Oncorhynchus masou) subspecies-complex (Salmonidae) is widely distributed in 
Far East Asia [1] and consists of four subspecies [2,3]. The subspecies are based on zoogeographic 
patterns and morphological characteristics; O. m. masou (Japanese name: sakura-masu), O. m. ishikawae 
(amago), O. m. subsp. (biwa-masu), and O. m. formosanus. O. m. masou and O. m. ishikawae have 
allopatric distributions, and are distinguishable by the spot-color pattern on their lateral sides: the former 
has black-spots while the latter has red-spots with sparse black-spots. O. m. subsp. is endemic to the 
Lake Biwa basin, central Honshu Island, Japan. Mitochondrial DNA analysis indicated that O. m. subsp. 
has a sister relationship with O. m. masou and O. m. ishikawae [4]. O. m. formosanus, which is endemic 
to Taiwan, has the world's southernmost distribution among the native salmonid fishes. 

The masu salmon subspecies-complex supports important recreational and commercial fisheries 
throughout the Japanese archipelago. Hatchery-reared masu salmon have been released into many rivers 
in Japan to enhance the fish stocks available to local fishermen and recreational anglers. However, their 
genetic integrities and diversities have been ignored, resulting in accelerated hybridizations between 
wild and hatchery fish within the subspecies [5,6] and among subspecies [7,8]. In addition, habitat 
fragmentations by artificial damming occur frequently in Japanese rivers, causing reduced population 
sizes and genetic diversities of salmonid populations [5,9,10]. O. m. subsp. is "near threatened (NT)" on 
the red list of the ministry of the environment, Japan (2007) [11]. O. m. formosanus is listed as an 
endangered species by the Taiwanese government, and as a critically endangered species by the 
International Union for Conservation of Nature (IUCN) [12-14]. 

Several dinucleotide-type microsatellite markers for masu salmon have been reported [15]. At 
present, however, we do not know whether these loci can be used across the masu salmon 
subspecies-complex. In addition, it will be useful to develop microsatellite markers with larger repeat 
motifs which allow easier geno typing than dinucleotide microsatellites. In this study, we performed 
454 pyrosequencing [GS-FLX+ System, 454 Life Sciences by Roche Diagnostics (Branford, CT, USA)] 
to isolate novel microsatellite markers from the O. m. masou genome, and evaluated their potential 
for conservation genetic applications, such as cross-subspecies amplifications and population 
genetic analyses. 

2. Results and Discussion 

Seven polymorphic microsatellite loci were isolated and deposited in GenBank 
(AB851460-AB851466; Table 1). The number of alleles detected at each locus ranged from 
four (OMAS-10) to 24 (OMAS-4) and the expected heterozygosity varied from 0.57 (OMAS-10) to 
0.92 (OMAS-5). The intrapopulation diversities differed, with the mean number of alleles per population 
varying between 4.43 (Miya River: O. m. ishikawae) and 11.29 (Chitose River: O. m. masou), and the 
expected heterozygosity varied between 0.60 (Miya River) and 0.80 (Chitose River; Table 2). Among 
27 Hardy- Weinberg equilibrium (HWE) tests, where OMAS-10 marker in O. m. ishikawae was omitted 
due to no allelic variation, six were significant after the false discovery rate correction of significant 
level for multiple comparisons. The observed departures, however, were neither locus nor population 
specific; therefore, we could not determine the cause of the HWE deviations. No marker pair showed 
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evidence of departure from linkage equilibrium, indicating no significant allelic association between the 
markers. These microsatellites represent useful tools for population genetic applications, such as stock 
identification in the Oncorhynchus masou complex. 



Table 1. Tetra-nucleotide microsatellite markers for Oncorhynchus masou complex. 



Marker 


Motif* 1 




Primer sequences (5'-3') * 2 


Ja* 3 


Ra (bp) * 4 


GenBank No. 


OMAS-3 


(AGAC)i4 


F-N 
R 


AGAGACAGATAGAGCCAGCCAG 
TGATGAACGTTACGATTGGAAG 


60 


127-199 


AB851460 


OMAS-4 


(CACT) 12 


F-P 
R 


TGCACATAAATTGAAGGCAAAC 
GAGTCACCTGTCCCTCAGTACC 


60 


137-249 


AB851461 


OMAS-5 


(TCTG) 13 


F-F 
R 


TTTTGCTGGTGACTCCCTAAAT 
ATTTGCAGAGGGAAACAGACAT 


60 


236-348 


AB851462 


OMAS-7 


(AGAC) 16 


F-N 
R 


GGGAAAGAAAGGAGATTGAGAGA 
GACCCTGGTAGAACTGCAAACT 


60 


235-295 


AB851463 


OMAS-10 


(GAGG) 6 


F-V 
R 


AGCAAAGGGAGATAAGGTAGGG 
CATCTTCATTCAGAGGGGTAGG 


60 


305-313 


AB851464 


OMAS-14 


(GAGT) 7 


F-F 
R 


ATTGTTAGAGCGGGAGACGATA 
TCCCCAGAATTGTTAGCTGAGT 


60 


98-138 


AB851465 


OMAS-18 


(ATCA) 9 


F-V 
R 


CACACAAAGAAAACCTTGAATGA 
AACGTGACTGCACAACAAAGTT 


60 


113-157 


AB851466 



* : Repeat counts are based on the individual from which the microsatellite-primer sequences were obtained; 
* 2 : The forward primer of each primer set was 5'-end-labeled with 6-FAM (denoted as F-F), VIC (F-V), 
NED (F-N) or PET (F-P) fluorescent dye; * 3 : Optimized annealing temperature. * 4 : Range of allele sizes. 



3. Experimental Section 

3.1. Isolation of Microsatellite Markers 

Total genomic DNA from O. m. masou (blood cell) was extracted by proteinase K digestion, 
phenol/chloroform extraction, and ethanol precipitation. We prepared a shotgun library from 1.0 \ig of 
genomic DNA following the 454-Roche protocol. Pyrosequencing was conducted on one-fourth 
picotiter plate. We searched for microsatellites and microsatellite-flanking PCR primer sequences 
using Auto-Primer [16]. A 454 run generated 175,689 sequence reads (63,019,585 bp) from ca. 
5.0 x 10 5 beads with clonally amplified DNA, resulting in 1937 contigs and 129,759 singletons. 
Auto-Primer proposed 242 PCR primer sets for tetra-nucleotide microsatellites. We selected 18 primer 
sets designed for amplification of "perfect-type" tetra-nucleotide repeats (i.e., repeat arrays without 
interruption by a non-repeat element); 1 1 were rejected because of poor amplification or difficulty in 
allele scoring. The remaining seven sets were applicable for estimation of polymorphism (Table 1). 
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Table 2. Genetic diversity statistics for the three subspecies of Oncorhynchus masou. 



Subspecies 


River/Lake 


Statistics of polymorphisms 1 


OMAS-3 


OMAS-4 


OMAS-5 


OMAS-7 


OMAS-10 


OMAS-14 


OMAS-18 






-"o 


0.86 


0.93 


0.73 


0.90 


0.50 


0.67 


0.31 








0.89 


0.90 


0.84 


0.86 


0.62 


0.75 


0.74 


Oncorhynchus masou masou 


Chitose River 


HW P 


0.12 


0.01 * 


0.00 * 


0.09 


0.10 


0.02 


0.00 * 






Mumnpr r\Y Allf^lf^o 
INUlilUtl Ul -rVllCICo 


i « 

i ~> 


1 7 


1 7 


1 ? 


j 


C 
o 


7 






n 


43 


42 


40 


42 


42 


43 


35 






0 


0.46 


0.83 


0.75 


0.24 




0.60 


0.15 






Hv 


0.44 


0.85 


0.75 


0.61 




0.44 


0.14 


Oncorhynchus masou ishikawae 


Miya River 


HW P 

i 1 VV 1 


1 00 


0 00 * 


0.02 


0 00 * 




0 0? 


1 00 

1 ,\J\J 






INUIIlUtl Ul -rVllClCo 


J 


Q 
o 




u 


1 


9 


"J 






n 


48 


48 


48 


41 


48 


48 


48 






H Q 


0.78 


0.87 


0.78 


0.78 


0.06 


0.53 


0.81 






H E 


0.81 


0.89 


0.86 


0.85 


0.06 


0.62 


0.78 


Oncorhynchus masou subsp. 


Lake Biwa 


HW P 

1 1 VV 1 


0 66 


0 80 


0 10 


0 1 1 


1 00 


0 60 


0 98 






Mumnpr r\T Allf^lf^o 
1NU111UC1 Ul -rlllClCo 


Q 


1 d 


1 1 


Q 


9 


7 


7 






n 


32 


31 


32 


32 


32 


32 


32 






H Q 


0.79 


0.64 


0.84 


0.71 


0.56 


0.79 


0.14 


Oncorhynchus masou masou 
(hatchery stock) 


Oohara River 


Hv 
ji t 

HW_P 

Number of Alleles 


0.71 
0.56 
7 


0.69 
0.47 
6 


0.86 
0.02 

8 


0.67 
0.69 
4 


0.53 
0.91 
3 


0.68 
0.37 
4 


0.27 
0.00 * 

3 






n 


39 


39 


38 


38 


39 


39 


37 


Average 




H 0 
He 


0.72 
0.71 


0.82 
0.83 


0.77 
0.83 


0.66 
0.75 


0.37 
0.40 


0.65 
0.62 


0.35 
0.48 






Number of Alleles 


9 


11 


11 


8 


2 


5 


5 



*', Ho: observed heterozygosity; H E : expected heterozygosity; HW P: probability value of Hardy- Weinberg equilibrium (HWE) testing. Significant HWE departure after the false discovery 
rate correction of significance level for multiple comparisons is indicated by an asterisk; n: number of genotyped individuals. 
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3.2. PCR Amplification and Genotyping 

Three wild population samples (O. m. masou from the Chitose River, Hokkaido Island, Japan; 
O. m. ishikawae from the Miya River, central Honshu, Japan; and O. m. subsp. from Lake Biwa, central 
Honshu, Japan) and one hatchery-reared stock sample (O. m. masou from the Ohara River, northern 
Honshu, Japan) were used for polymorphism detection. After conjugating 6-FAM, HEX, NED or PET 
fluorescent dyes (Applied Biosystems, Life Technologies, Carlsbad, CA, USA) to the 5'-end of each 
forward primer, we performed multiplex PCR for the seven selected primer sets using a Qiagen 
multiplex PCR kit (Qiagen, Limburg, The Netherlands). The 10-|xL reaction mixture contained 
1 x Qiagen multiplex PCR master mix, 0.2 \iM of each primer, 2 jxL of distilled water, and 1 |xL of DNA 
solution. Amplifications were carried out in a GeneAmp PCR System 2700 thermal cycler (Applied 
Biosystems), according to the supplier's instructions (Qiagen multiplex PCR kit; Qiagen): initial 
denaturation at 95 °C for 15 min; followed by 30 cycles of denaturation at 94 °C (30 s), annealing (60 °C, 
90 s) and extension (72 °C, 60 s); with a final extension (60 °C, 30 min). PCR fragments and GeneScan 
600LIZ size standards were resolved on an ABI PRISM 3730x/ sequencer (Applied Biosystems, 
Life Technologies). Allele size was determined using GeneMapper v4.1 (Applied Biosystems, 
Life Technologies). 

3.3. Data Analysis 

Polymorphism statistics, including the observed/expected heterozygosities, were calculated in 
Arlequin v3.5 [17]. The Hardy- Weinberg equilibrium of markers/populations was tested using Fisher's 
exact test with the Markov chain method (Markov chain steps: 10 5 ; Dememorization: 10 5 ). We also 
examined linkage equilibrium between markers (likelihood-ratio test: 10 4 permutations). Critical 
significance levels for multiple testing were adjusted following the false discovery rate [18]. 

4. Conclusions 

We report the isolation and characterization of seven polymorphic microsatellite loci for the 
Oncorhynchus masou complex. All seven markers used were polymorphic. Cross-subspecies 
amplification for O. m. masou, O. m. ishikawae and O. m. subsp. was successful. These microsatellites can 
be utilized in studies of genetic structure, genetic diversity, and intra- and inter-subspecific hybridization, 
making a contribution to conservation and management of the Oncorhynchus masou complex. 
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